A new class of exponential propagation iterative methods of Runge-Kutta type (EPIRK)
نویسنده
چکیده
We propose a new class of the exponential propagation iterative methods of RungeKutta-type (EPIRK). The EPIRK schemes are exponential integrators that can be competitive with explicit and implicit methods for integration of large stiff systems of ODEs. Introducing the new, more general then previously proposed, ansatz for EPIRK schemes allows for more flexibility in deriving computationally efficient highorder integrators. Recent extension of the theory of B-series to exponential integrators [1] is used to derive classical order conditions for schemes up to order five. An algorithm to systematically solve these conditions is presented and several new fifth order schemes are constructed. Several numerical examples are used to verify the order of the methods is verified and illustrate the performance of the new schemes.
منابع مشابه
New Adaptive Exponential Propagation Iterative Methods of Runge-Kutta Type
Exponential integrators have emerged as an efficient alternative to commonly used time-integrators. Recently a new class of exponential propagation iterative methods of Runge-Kutta type (EPIRK) has been introduced [30]. These schemes possess a structure that makes them computationally advantageous compared to other exponential methods. In addition, the general EPIRK formulation offers flexibili...
متن کاملEPIRK-W and EPIRK-K time discretization methods
Exponential integrators are special time discretization methods where the traditional linear system solves used by implicit schemes are replaced with computing the action of matrix exponential-like functions on a vector. A very general formulation of exponential integrators is offered by the Exponential Propagation Iterative methods of Runge-Kutta type (EPIRK) family of schemes. The use of Jaco...
متن کاملSolving Parameter Estimation Problems with Discrete Adjoint Exponential Integrators
The solution of inverse problems in a variational setting finds best estimates of the model parameters by minimizing a cost function that penalizes the mismatch between model outputs and observations. The gradients required by the numerical optimization process are computed using adjoint models. Exponential integrators are a promising family of time discretizations for evolutionary partial diff...
متن کاملExponential Runge-Kutta Methods for Stiff Kinetic Equations
Abstract. We introduce a class of exponential Runge-Kutta integration methods for kinetic equations. The methods are based on a decomposition of the collision operator into an equilibrium and a non equilibrium part and are exact for relaxation operators of BGK type. For Boltzmann type kinetic equations they work uniformly for a wide range of relaxation times and avoid the solution of nonlinear ...
متن کاملRunge-Kutta convolution quadrature for operators arising in wave propagation
An error analysis of Runge-Kutta convolution quadrature is presented for a class of nonsectorial operators whose Laplace transform satisfies, besides the standard assumptions of analyticity in a half-plane Re s > σ0 and a polynomial bound O(s 1) there, the stronger polynomial bound O(s2) in convex sectors of the form | arg s| ≤ π/2 − θ < π/2 for θ > 0. The order of convergence of the Runge-Kutt...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- J. Comput. Physics
دوره 230 شماره
صفحات -
تاریخ انتشار 2011